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Resumen 


Tzatchkov, V. G., € Alcocer-Yamanaka, V. H. (mayo-junio, 
2016). Modelación de la variación del consumo de agua 
potable con métodos estocásticos. Tecnología y Ciencias del 
Agua, 7(3), 115-133. 


A pesar de su importancia, en la práctica, la variación de 
la demanda de agua potable se estima de manera muy 
aproximada de mediciones continuas del gasto en la tubería 
que abastece a una zona; por su lado, la curva de variación 
de la demanda medida se asume válida para cualquier 
otra tubería de la red, independientemente del número de 
usuarios al que da servicio. El presente artículo describe una 
metodología para obtener la variación diaria estocástica de 
la demanda instantánea de agua potable, aplicable a una sola 
casa o cualquier número de casas con base en los parámetros 
estadísticos del consumo de agua en casas individuales, 
nivel de fugas dado y la variación del gasto medida en la 
tubería de abastecimiento, considerando casos de suministro 
continuo e intermitente. La metodología propuesta ofrece 
una base racional para la determinación de la variación de 
la demanda a cualquier nivel de agregación, que puede 
combinarse con la correspondiente manera tradicional o 
sustituirla. Los resultados se compararon con mediciones de 
campo en una ciudad mexicana donde el suministro de agua 
potable es continuo. La variación de los gastos en las tuberías 
con suministro intermitente o con servicio continuo, pero 
con cisternas y tinacos en los domicilios, es muy diferente 
de aquella en tuberías con suministro continuo sin cisternas 
y tinacos. El coeficiente de demanda máxima horaria es más 
alto en suministro intermitente que en suministro continuo, 
pero más bajo en redes que tienen servicio continuo, y 
cisternas o tinacos en los domicilios. 


Palabras clave: demanda estocástica de agua potable, 
método de pulsos rectangulares de Poisson (PRP), curva de 
la variación horaria de la demanda, coeficiente de variación 
horaria, agregación de la demanda de agua, suministro 
continuo e intermitente de agua potable. 


Abstract 


Tzatchkov, V. G., € Alcocer-Yamanaka, V. H. (May-June, 2016). 
Stochastic Method Water Demand Variation Modelling. Water 
Technology and Sciences (in Spanish), 7(3), 115-133. 


Despite its importance, water demand variation in current practice 
is estimated very approximately by observing water flow in a pipe 
supplying certain zone, and the demand variation curve thus 
obtained is assumed to be valid for any other pipe in the network, 
regardless the number of water users served. This paper describes 
a methodology for obtaining the daily instantaneous water demand 
stochastic variation based on statistical parameters of water 
consumption at individual homes, given leakage level, and observed 
flow variation at the supplying pipe. Continuous and intermittent 
supplies are considered. The methodology is applicable to a single 
home or a group of any number of homes. It provides a rational base 
for determining water demand variation at any aggregation level 
and can be combined with the corresponding traditional way of such 
determination, or substitute it. The results were compared with field 
measurements in a Mexican city where water supply is continuous. 
Flow variation in intermittent water supply pipes, or in continuous 
supply pipes with ground level cisterns and roof tanks, is very 
different from that in continuous supply pipes. The hourly water 
demand coefficient is higher in intermittent water supply, compared 
to continuous supply, but lower in networks with continuous supply 
and ground level cisterns and roof tanks. 


Keywords: Stochastic water demand, PRP method, Hourly demand 
variation curve, Hourly demand variation coefficient, Water demand 
aggregation, Continuous and intermittent water supply. 
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Introducción 


La estimación correcta de la demanda de agua 
potable representa una condición indispensable 
para la planeación y el diseño de los sistemas de 
suministro, que en gran medida determina las 
inversiones necesarias y calidad del servicio. La 
demanda está sujeta a variaciones interanuales, 
estacionales, semanales y diarias. La variación 
interanual se expresa por la evolución de la 
demanda durante varios años debida a la 
dinámica de la población, alteraciones en su 
estilo de vida y, recientemente, por el cambio 
climático. La variación estacional se influye por 
las diferencias en el clima en las estaciones del 
año; y la variación semanal muestra carácter 
cíclico con cierto patrón diario, que es similar 
entre un día y otro, con algunas diferencias 
para los fines de semana y días festivos. A pesar 
de su importancia, en la práctica, la demanda 
de agua potable y su variación se estima de 
manera muy aproximada. En la ingeniería 
del abastecimiento de agua potable es común 
utilizar los conceptos coeficientes de variación 
diaria y horaria de la demanda y curva de variación 
horaria de la demanda. Llamado en algunos países 
también coeficiente pico o coeficiente punta, el 


coeficiente de variación de la demanda expresa, 
para una tubería dada, la relación entre el gasto 
máximo y el gasto medio que conduciría la 
tubería, que se utiliza ante todo para revisar la 
capacidad de la tubería en el diseño. La curva 
de variación horaria de la demanda para una lo- 
calidad normalmente se obtiene de mediciones 
continuas del gasto en una tubería que sale de 
una fuente o un tanque para abastecer la red de 
distribución. Una vez obtenidos de esta manera 
para cierta tubería, los coeficientes de variación 
y curva de variación de la demanda se asumen 
válidos para cualquier otra tubería de la red, con 
independencia del número de usuarios a que da 
servicio, e incluso para los propios domicilios. 
La figura 1 muestra la curva de variación de la 
demanda que maneja la Comisión Nacional del 
Agua (MAPAS-Conagua 2007a), obtenida para 
varias poblaciones de México. 

En el ámbito domiciliario, la demanda real 
de agua potable —definida por las necesidades 
y costumbres de sus habitantes—dista mucho, 
sin embargo, de la representación dada por la 
figura 1. La figura 2 muestra el consumo de 
agua registrado en un día en una casa mexica- 
na con resolución ultra fina (cada segundo). El 
consumo se presenta con pulsos instantáneos, 
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Figura 1. Curva de la variación horaria de la demanda de agua potable para la república mexicana, recomendada en los 


manuales de la Comisión Nacional del Agua (MAPAS-Conagua, 2007a). 
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por lo general de poca duración, seguidos por 
periodos prolongados sin consumo. En un 
estudio previo (Izatchkov, Alcocer, Arreguín, 
éz Feliciano, 2005; Alcocer, 2007), se obtuvo 
que sólo de 3 a 5% del total del día se destina 
a actividades de consumo de agua. Dado que 
los momentos de tiempo exactos en que los 
usuarios hacen uso del agua no serán los mis- 
mos entre un día y otro, y entre una casa y otra 
casa, la variación del consumo tiene también 
un carácter estocástico, sobre todo en escalas de 
tiempo finas del orden de segundos. El consumo 
total de un pequeño grupo de domicilios sería 
similar al mostrado en la figura 2, y tenderá al 
mostrado en la figura 1 para un gran número 
de usuarios en que también se incluyen fugas. 

Aparte de la determinación del coeficiente 
de variación de la demanda, que se emplea para 
revisar la capacidad de las tuberías en el diseño, 
es necesario considerar la variación temporal de 
la demanda de agua potable en al menos otros 
dos tipos de análisis: 


a) En el análisis del balance de las masas de 
agua que entran y salen de un tanque, para 
determinar o revisar el volumen de almace- 
namiento necesario. 


b) En una simulación (modelación) dinámica 
del comportamiento hidráulico y de calidad 
de agua en el sistema de distribución (esta 
simulación puede incluir también el punto 
(a) anterior). En las últimas décadas se 
han realizado avances importantes en la 
modelación matemática de redes de dis- 
tribución de agua. Existen ahora sistemas 
computacionales (software) de libre acceso, 
como Epanet (Rossman, 2002), o comerciales, 
como InfoWorksG (Innovyze, 2015), y otros 
(MAPAS-Conagua, 2007b), que permiten 
introducir la red completa (incluyendo red 
secundaria) en un ambiente gráfico. Algu- 
nos de estos sistemas pueden integrarse 
con sistemas de información geográfica, 
con grado de detalle que puede incluir 
hasta cada toma domiciliaria en el modelo 
(Rodríguez et al., 2010). Con independencia 
del grado de detalle del modelo (un nodo 
por toma o un nodo por grupo de tomas), 
el modelo debe ser alimentado con la varia- 
ción temporal de la demanda en cada nodo, 
que a su vez depende de la cantidad de 
tomas servidas por el nodo y es estocástica. 
En los modelos de redes de agua potable 
actuales esto normalmente no se hace, 
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Figura 2. Variación típica del consumo de agua medido en un día (Alcocer, 2007). 
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asumiendo una curva de variación horaria 
igual y determinista para todos los nodos. 
La modelación de las redes con el uso de 
ese tipo de curvas resulta aceptable para 
las tuberías principales, donde la variación 
del caudal conducido es continua, pero es 
poco realista en las tuberías secundarias que 
proporcionan el servicio, generando sobre 
o sub dimensionamiento en los diseños, e 
incertidumbre en los resultados de los mo- 
delos de redes existentes, que se expresa en 
discrepancias entre los valores medidos y 
calculados de gastos y presiones. 


Ante esta situación, para modelar con 
mayor certeza el comportamiento hidráulico 
y de calidad del agua dentro de las redes de 
distribución, incluyendo tuberías secundarias, 
surge la necesidad de enfoques más realistas, 
que incluyan una definición temporal precisa 
del consumo doméstico con escalas inferiores a 
la horaria, considerando su carácter estocástico. 
Para responder a esta necesidad, en el presen- 
te artículo se propone una metodología que 
permite obtener la variación diaria estocástica 
de la demanda instantánea de agua potable a 
diferentes niveles de agregación temporal, apli- 
cable a una sola casa o un grupo de N casas. Con 
el uso de esta metodología se hace un nuevo 
planteamiento del concepto y la obtención del 
coeficiente de variación de la demanda, que 
ahora se considera dependiente de la cantidad 
de usuarios abastecidos. 


Metodología 


Existen técnicas con enfoque estocástico para 
representar el consumo doméstico de agua, 
desarrolladas básicamente en los últimos años, 
tales como el Método de Pulsos Rectangulares 
de Poisson (en inglés, Poisson Rectangular Pulse 
Method, abreviado como PRP method) (Buch- 
berger £ Wu, 1995; Buchberger, Carter, Lee, e 
Schade, 2003) y el Método de Pulsos Rectan- 
gulares de Neyman-Scott (en inglés, Neyman- 
Scott Rectangular Pulse Method, abreviado 
como NSRPM; Neyman éz Scott, 1958.) (Alvisi, 


+» ISSN 0187-8336 


Franchini, €: Marinelli, 2003; Alcocer ef al., 2008). 
Estos dos métodos manejan la demanda de agua 
y los parámetros para modelarla, en el nivel 
toma domiciliaria, sin necesitar información de 
los aparatos sanitarios que se tienen dentro de 
las casas. El método PRP maneja directamente 
la demanda instantánea (con intervalo de me- 
dición, registro y modelación de un segundo), 
mientras que NSRPM puede trabajar con dife- 
rentes intervalos de registro, al considerar una 
desagregación temporal de la demanda (Alcocer 
et al., 2008). En los trabajos de Blokker (Blokker, 
2010; Blokker, Pieterse-Quirijns, Vreeburg, € 
Van Dijk, 2011), el método PRP se extiende a 
modelar, empleando el método de Monte Carlo, 
los pulsos de demanda que generan los aparatos 
sanitarios individuales en cada domicilio, que 
luego son agregados al nivel deseado para 
obtener los gastos en la tubería que suminis- 
tra a un grupo de domicilios. Este método es 
conceptualmente más completo y exacto, pero 
requiere de datos detallados de la cantidad, tipo, 
frecuencia, patrón de uso y capacidad de todos 
los aparatos sanitarios en cada domicilio. En 
una comparación entre su método, que es im- 
plementado en un software llamado SIMDEUM, 
Blokker, Buchberger, Vreeburg y Van Dijk (2008) 
reportan que se necesitaron 99 parámetros para 
modelar la demanda con SIMDEUM, compara- 
do con sólo cinco parámetros en el método PRP. 
Esto lo vuelve impráctico, en particular para 
países como México, donde es problemático 
obtener información de lo que hay dentro de 
los domicilios. 

Magini, Pallavicini y Guercio (2008), y Ver- 
tommen, Magini y Cunha (2012) enfocan el pro- 
blema del carácter estocástico y la incertidumbre 
de la demanda de agua potable por medio de 
escalamiento temporal y especial. Otros autores 
que han tratado la modelación de la demanda 
estocástica y sus efectos son García, García, 
Cabrera, Arregui y García (2004) (que usan el 
método PRP), y Van Zyl (Van Zyl, 1996; Van Zyl, 
le Gat, Piller, £: Walski, 2012), cuyo trabajo ha 
sido dedicado ante todo a los efectos sobre la 
confiabilidad de los sistemas de agua potable. 
Alcocer, Tzatchkov y Arreguín (2012) propu- 
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sieron una metodología para la modelación de 
redes de agua potable con demanda estocástica 
y compararon los resultados con mediciones en 
una red real. 


Desarrollo del presente estudio 


Para lograr el propósito planteado, el estudio se 
desarrolló de la manera siguiente: 


1. Se analizaron los métodos existentes de 
modelación de demanda estocástica de 
agua potable, considerando la información 
necesaria para su implementación, y la 
eventual disponibilidad de esta información 
O facilidad para obtenerla. Se escogió en 
resultado el Método de Pulsos Rectangu- 
lares de Poisson (PRP), que ocupa menos 
parámetros y requiere de menos datos. 

2. Se definieron los parámetros que caracteri- 
zan la variación del consumo en el método 
PRP para los fines del estudio, y las ecua- 
ciones que los relacionan, considerando su 
carácter estocástico. 

3. Se desarrolló un procedimiento de simu- 
lación numérica basado en este método, 
que se implementó en Visual Basic (VBA) 
para Excel, proporcionando la variación 
estocástica de la demanda de agua potable 
a diferentes niveles de agregación y los 
coeficientes de variación en función de la 
población servida. 

4. El procedimiento propuesto se aplicó a una 
ciudad mexicana donde el suministro de 
agua potable es continuo. Los resultados se 
compararon con mediciones de campo y se 
obtuvieron los coeficientes de variación de 
la demanda. 

5. El análisis de los coeficientes de variación se 
amplió a diferentes niveles de agregación, y 
casos de suministro continuo y suministro 
intermitente, tomando en cuenta las fugas. 
Se consideró tan sólo el ciclo diario y los 
componentes aleatorios, dejando de lado 
los componentes interanuales, estacionales 
y semanales, que no son puramente estocás- 
ticos y pueden ser considerados aparte. 


Parámetros que caracterizan la variación del 
consumo 


El modelo PRP asume que la demanda se 
presenta en pulsos rectangulares (en el plano 
tiempo-demanda), con duración e intensidad 
aleatorias, cuyo surgimiento en el tiempo obe- 
dece a un proceso de Poisson no homogéneo. 
La aplicación de ese modelo requiere de varios 
parámetros básicos de la demanda local en los 
domicilios, como la frecuencia de uso del agua 
por hora del día (u otro intervalo de tiempo), 
intensidad y duración promedio de los pulsos, 
y sus coeficientes de variación. Buchberger et 
al. (2003), Alcocer, Tzatchkowv Buchberger y 
Feliciano (2004), y Tzatchkov et al. (2005) obtu- 
vieron parámetros de este tipo para la demanda 
de agua potable en localidades de los Estados 
Unidos de América y México. 

El proceso de Poisson describe las series de 
ocurrencias independientes en el eje del tiempo 
mediante el mecanismo estocástico, separadas 
entre sí por intervalos cuya distribución de 
probabilidades es exponencial, dada por la 
ecuación: 


f(t)=h2e" (1) 


donde A expresa la frecuencia media con que 
ocurren los pulsos de demanda; Á se mide en 
pulsos por unidad de tiempo, y el contexto del 
método PRR se conoce por lo común con el 
nombre tasa de llegada. 

La función de distribución cumulativa co- 
rrespondiente a la distribución (1) es: 


F(t)=1-e"" (2) 
donde por definición O < F(t) < 1. El valor medio 
es igual a 1/4. Señalando F(t) con x, la solución 
de la ecuación (2) con respecto al tiempo t es: 


1 
¿== In(1-x) (3) 


El proceso es aditivo (Buchberger et al., 
2003), por lo que para un grupo de N casas de 
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parámetros de demanda de agua similares, la 
ecuación (3) se convierte en: 


La ecuación (5) proporciona el tiempo entre 
las ocurrencias (pulsos de demanda de agua), 


correspondiente a cierta probabilidad x, y 


e - in (1 Ñ x) (4) permite obtener este tiempo en el proceso de si- 


mulación para valores de x (x > 0) generados de 


; s : : una distribución de probabilidades uniforme. La 
Si x es un número aleatorio obtenido de la ñ 3 ' na 
distribución uniforme U(0,1), entonces 1 — x gula A PA 
a la ecuación (5) para N.A=3 600 pulsos / hora (un 
también pertenece a U(0,1), con lo que la ecua- Ñ j 
vá e ; pulso /segundo), que podría ser obtenido por 
ción (4) se puede simplificar como: 

3 600 casas con un pulso por hora en cada casa, u 


1 otra combinación entre N.A. La figura 4 muestra 


t= a (5) lo mismo para para N.A =4 pulsos / hora. 
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Figura 3. Tiempo entre pulsos de demanda para diferentes valores del número aleatorio, para tasa de llegada 
de un pulso por segundo. 
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Figura 4. Tiempo entre pulsos de demanda para diferentes valores del número aleatorio, para tasa de llegada 
de cuatro pulsos por hora. 
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Procedimiento de simulación del consumo de 
agua potable 


Se desarrolló un procedimiento de simulación 
del consumo de agua potable en un domicilio, o 
en un grupo de domicilios, basado en el método 
de los pulsos rectangulares de Poisson (Método 
PRP). Se considera que los pulsos de demanda 
ocurren en instantes aleatorios, y cada pulso 
tiene una duración e intensidad aleatorias. La 
tasa de llegada varía durante el día, por lo que 
se tiene un proceso de Poisson no homogéneo. 
Los datos que se requieren para la simulación 
son los siguientes: 


e Intervalo de tiempo que se quiere modelar; 
por lo normal se maneja un intervalo de un 
día. 

+ Tasa de llegada de los pulsos (pulsos / hora) 
por cada hora (u otro intervalo dado) del 
día. 

+ Duración media de los pulsos (s). 

e Distribución de probabilidades de la 
duración de los pulsos. Las investigaciones 
realizadas por Buchberger et al. (2003) 
mostraron que la duración de los pulsos 
se describe de forma satisfactoria por una 
distribución log normal, por lo que se usó 
ese tipo de distribución probabilística. 

e Desviación estándar de la duración de 
los pulsos (s). Dado que se emplea una 
distribución log normal, debe ser expresada 
por el logaritmo natural (In x) de dicha 
duración (s). 

+ Intensidad media de los pulsos (1/s). 

e Distribución de probabilidades de la 
intensidad de los pulsos. Las investigaciones 
realizadas por Buchberger et al. (2003) 
mostraron que la intensidad de los pulsos 
se describe también por una distribución 
log normal, por lo que se usó ese tipo de 
distribución probabilística. 

e Desviación estándar de la intensidad de los 
pulsos (1/s). De manera similar, dado que se 
emplea una distribución log normal, debe 
ser expresada por el logaritmo natural (In 
x) de dicha intensidad (1/s). 


e Cantidad de casas. 


En el caso de un proceso de Poisson no 
homogéneo (es decir, con A variable en el 
tiempo), la aplicación directa de la ecuación (5) 
produciría resultados incorrectos, dado que A 
ya tendría otro valor para el momento t que se 
busca; es decir, se despreciaría la variación de 
A en el intervalo de tiempo entre los dos pulsos 
consecutivos. Se han propuesto varios métodos 
para tratar este problema (Leemis € Park, 2006), 
entre los cuales se seleccionó el de inversión de 
la función de la tasa de llegada acumulada, 
planteado por Cinlar (1975). Sea la tasa de 
llegada de los pulsos de demanda 1, variable 
durante el periodo T a modelar, representada 
por medio de k+ 1 puntos (t, A,), con O = t,<t, 
<... <t, =T, donde h, > 0. Se asume que A varía 
linealmente entre cada punto t, y t,.,. La tasa de 
llegada acumulada se define como: 


A¡= f/A(t)dt — j=0,L...k (6) 


y se calcula recursivamente como: 
A, =0 


Atl) 0 


La pendiente de cada segmento s, se calcula 
como: 

MA, 
ia A ¡=0,1,...,k (8) 


i 
bj Ey 


Para el subintervalo tstst,, 


A+) =2,+s,(+-,) (9) 


M1) A, +2, (t=4)+35/(1-4,) (10) 


Ahora bien, para obtener los momentos en 
que se producen los pulsos de demanda, se 
genera un proceso auxiliar de Poisson homogé- 
neo, en el intervalo de 0 a 1, con tasa de llegada 
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unitaria, es decir, con A = 1 en la ecuación (3). 
Este proceso obtendrá una serie de sucesos u,, 
Uy... U,. Luego, para cada u se obtendrá el ins- 
tante en que se produce un pulso de demanda 
despejando t de la ecuación (10): 


(11) 


Is Au-A) 


á A+ pi +2s,(u—A,) 


La figura 5 ilustra este proceso. 
En términos algorítmicos, este proceso se 
describe de la manera siguiente: 


1. Mediante integración numérica (con regla 
trapezoidal) se obtiene la función que re- 
presenta la tasa de llegada acumulada de 
los pulsos para cada hora (u otro intervalo 
dado) del día. 

2. Inicializar t,= 0.0; u,=0.0; n =0. 

3. Mientras f,< TT. 

a) Generar un número aleatorio x entre 0 
y 1 basado en distribución de probabi- 
lidades uniforme. 

b) Calcular con este número aleatorio, por 
la ecuación (5) y A = 1, el lapso en que 
ocurre el siguiente suceso, es decir, y, 
=4,-Inx. 

c) Obtener el instante en que inicia el pulso 
de demanda siguiente de la ecuación 
(11). 


1 


d) Generar otro número aleatorio x entre 0 
y 1 basado en distribución uniforme. 

e) Evaluar la función inversa de la distri- 
bución log normal con este número. El 
resultado, un número aleatorio de dis- 
tribución de probabilidades log normal, 
se multiplica por la duración media de 
los pulsos, obteniéndose la duración 
aleatoria del pulso. Se ejecuta el mismo 
proceso para lograr la intensidad alea- 
toria del pulso. 

f) Se registra el pulso con su instante de 
inicio (en segundos), duración e inten- 
sidad así conseguidos; es decir, se suma 
la demanda con la intensidad obtenida 
en los tiempos que corresponden a su 
momento de inicio f, y su duración. 

Fin del ciclo. 
Como resultado de los puntos 1 a 4 del 
algoritmo se obtiene una serie de demanda 
estocástica para el intervalo deseado que 
tiene resolución temporal de un segundo. 
Se calculan los parámetros estadísticos de 
esta serie, como media y varianza. 
Para comparar el resultado obtenido con 
series medidas, que normalmente no tienen 
resolución de un segundo, se obtienen los 
promedios de la serie para intervalos desea- 
dos, por ejemplo, para cada minuto o cada 
hora. 


Uy 


U3 


u 


Figura 5. Obtención de los tiempos en que se producen los pulsos de demanda por medio de un proceso de 


Poisson homogéneo auxiliar. 
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Los datos que se ocupan en el procedimiento 
propuesto, es decir, los parámetros estadísticos 
de la duración e intensidad de los pulsos de 
demanda, así como su tasa de llegada por cada 
hora del día, son en principio diferentes para 
cada lugar y pueden conseguirse por medición 
directa en cierto número de tomas que abastecen 
a domicilios. La metodología correspondiente 
de esta medición está descrita en Alcocer et al. 
(2004), Tzatchkov et al. (2005), y Alcocer (2007), 
con resultados para una ciudad mexicana. Estas 
mediciones se pueden complementar con medi- 
ciones del gasto suministrado en las fuentes de 
abastecimiento para determinar el porcentaje de 
fugas de la siguiente manera: 


+ El producto de la duración media e intensi- 
dad media proporciona el volumen de agua 
promedio de un pulso de demanda. 

+ El volumen promedio del pulso de deman- 
da se multiplica por la cantidad de pulsos 
en una hora para obtener el volumen con- 
sumido en los pulsos de demanda durante 
la misma hora. 

e La diferencia entre el volumen de agua re- 
gistrado en una hora y el volumen consumi- 
do proporciona el volumen de fugas. 


De otro modo, si se conoce el porcentaje de 
fugas, la medición del gasto suministrado en las 
fuentes se puede usar para ajustar el número 
de pulsos por hora (la tasa de llegada A) de la 
siguiente manera: 


+ Del volumen de agua registrado en una 
hora se restan las fugas. El resultado es 
el volumen consumido en los pulsos de 
demanda durante la misma hora. 

e El volumen consumido en los pulsos de de- 
manda se divide entre el volumen de agua 
promedio de un pulso de demanda para ob- 
tener la cantidad de pulsos en la hora. 


Aplicación a un sistema real 


A manera de ejemplo, se presenta la aplicación 
del procedimiento propuesto a una zona aislada 
de Tijuana, Baja California, abastecida de una 
sola fuente, para la cual se tienen mediciones 
del gasto en la fuente, con resolución de 60 s. La 
figura 6 muestra la variación medida del gasto 
en un día típico. El servicio de agua en la zona es 
continuo y las casas no tienen tinacos o cisternas. 
El gasto medio medido durante este día es de 
33.01 1/s, correspondiente a unas 3 300 casas. 
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Figura 6. Variación medida del gasto en la tubería que abastece una zona aislada. 
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Los parámetros estadísticos de los pulsos 
de demanda, tomados de un estudio anterior 
de caracterización de la demanda estocástica 
en domicilios mexicanos con servicio continuo 
(Alcocer et al., 2004; Alcocer, 2007; Tzatchkov et 
al., 2005) se muestran en el cuadro 1. 

La tasa de llegada de los pulsos, determinada 
con el procedimiento arriba descrito, asumiendo 
un 18% de fugas, se representa gráficamente en 
la figura 7. 

La figura 8 muestra la serie de demanda 
generada por el procedimiento propuesto con 
anterioridad para la misma cantidad de 3 300 
casas. Esta serie tiene una resolución temporal 
de un segundo. Tratándose de un proceso esto- 
cástico, con cada corrida del procedimiento se 
obtendrá una serie de demanda diferente, que 
en los términos de la modelación estocástica se 
llamaría una realización del proceso. Las reali- 


zaciones obtenidas deben exhibir los mismos 
parámetros estadísticos de las series generadas, 
tales como media y desviación estándar. 

Una vez generada una serie sintética, la 
misma puede ser agregada a cualquier inter- 
valo deseado. La figura 9 muestra la serie de 
demanda de la figura 8, agregada a 60 segundos, 
y comparada con el gasto medido, ambas en 
forma adimensional. Tratándose de un proceso 
estocástico, la coincidencia se puede calificar de 
muy buena a excelente. La figura 10 muestra la 
serie agregada a una hora, que replica aproxi- 
madamente la variación dada de la tasa de 
llegada de los pulsos. 

Una ventaja importante del procedimiento 
desarrollado es que permite generar series para 
diferente número de casas. Una vez obtenida la 
serie de demanda para la cantidad de casas para 
la cual se dispone de medición de la variación 


Cuadro 1. Parámetros de los pulsos de demanda. 


Duración media de los pulsos (s) 74.4 
Desviación estándar de LN (duración) (s) 0.46 
Intensidad media de los pulsos (1/s) 0.1453 
Desviación estándar de LN (intensidad) (1/s) 0.1855 
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Figura 7. Tasa de llegada de los pulsos de demanda. 
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Figura 9. Serie de demanda generada en forma adimensional, agregada a 60 segundos, comparada con el gasto medido. 


del gasto, y verificada por comparación con la 
serie medida, con el procedimiento propuesto 
se pueden obtener series de demanda para 
cualquier cantidad de casas, para su uso en 


modelos dinámicos de redes de distribución o 
para otros fines. La figura 11 muestra la serie 
generada para una sola casa, la figura 12 para 20 
casas, y la figura 13 para 100 000 casas. La serie 


Tecnología y Ciencias del Agua, vol. VIL, núm. 3, mayo-junio de 2016, pp. 115-133 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VU, núm. 3, mayo-junio de 2016, pp. 115-133 


Tzatchkov £ Alcocer-Yamanaka, Modelación de la variación del consumo de agua potable con métodos estocásticos 


Serie agregada a 3 600 s 
60 
50 
_w 
= 
3 
E 30 
a] 
E 
A 
20 
10 
0 e e e e e e e e e e e e 
[a] Q [=] [a] e [a] Q e [a] e o Q 
N ÁH o) [o] o NA + No 59) [a] N X= 
[53 + Sl o No cm o 133 $ AN key] No 
+ A AN cm il LO LO No Ñn 153 00 
Tiempo (s) 
Figura 10. Serie de demanda generada, agregada a una hora. 
0.35 
0.3 
0.25 
Edo 
= 0.2 
a 
ll 
E ] 
E 0.15 
A 
085 =] = e Sa = =] e = e =] 2 e 
e Q [==] [e] e e Q [e] [e] e [e] Q 
SS + No] ES S SS + No) ES S A + 
ES = ES Si 3 4 8 ÍS NS IS ES 
Tiempo (s) 


Figura 11. Serie de demanda generada para una casa. 


para 100 000 casas, sin ser una serie agregada, 
prácticamente replica la forma de la curva que 
representa la tasa de llegada de los pulsos, que 
es de esperarse, pues para gran número de casas 
(para N grande), la ecuación (4) da tiempos en- 
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tre la llegada de pulsos incluso mucho menores 
de un segundo. 

A pesar del procedimiento estocástico de 
generación de los pulsos de demanda, no se pro- 
ducen valores excesivos de su intensidad para 


Tzatchkov € Alcocer-Yamanaka, Modelación de la variación del consumo de agua potable con métodos estocásticos 


1.2 


o 
00 R 
on 


Demanda (1/s) 
o 
R 


IN 
ON "5 ] Ú 
' A IM 


e 
a 


0.2 


MN mi 


28 800 
36 000 
43 200 === 
50 4004 
57 600 pra 


A e e e e 
[a] S e [a] [a] 
+ 00 o A Ál 
= 3 OR RS 
Tiempo (s) 
Figura 12. Serie de demanda generada para 20 casas. 
1800 
1600 
1.400 
1200 
Ñ 
= 1000 
ls] 
Mo] 
E 800 
E 
A 
600 
400 
200 
A 
e Q [a] [e] e [e] Q [»] e e e Q 
A + No) 00 o A + No) 00 o A + 
3 + a! o No] ca o 55 y AN DD Noel 
al AN AN a al LO [Lo] NO Ñn 153 o 
Tiempo (s) 
Figura 13. Serie de demanda generada para 100 000 casas. 
un número pequeño de domicilios, como puede 2007), que pueden requerir de restricciones 
suceder en otros métodos de simulación de la adicionales para limitar tales valores excesivos. 
demanda estocástica, por ejemplo en el NSRPM Esto es otra ventaja del método propuesto en 
(Alvisi et al., 2003; Alcocer et al., 2008; Alcocer, este artículo. 
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Coeficientes de variación de la demanda 


El método de modelación de las series de de- 
manda propuesto hace posible obtener también 
estimados para el coeficiente de variación de 
la demanda, dado por la relación del gasto 
máximo producido entre el gasto medio (no 
confundir con el coeficiente de variación de la 
serie que se maneja en estadística, que se define 
como la relación entre desviación estándar y 
media aritmética). En la práctica de la ingeniería 
hidráulica se maneja el coeficiente de variación 
horaria de la demanda, pero el método de mode- 
lación de las series de demanda propuesto hace 
posible que se obtenga para cualquier intervalo 
deseado. Este coeficiente será tanto menor cuan- 
to mayor sea el periodo de agregación temporal, 
debido al promediado en la agregación, y tanto 
menor cuanto mayor sea la cantidad de casas a 
que corresponde la serie de demanda. La figura 


14 compara la serie medida de gasto cada 60 
s del ejemplo anteriormente presentado (con 
gasto medio de 33.01 1/s, correspondiente a 
3 300 casas) y su promediado horario. La figura 
15 realiza la misma comparación para una serie 
de gasto medido cada 60 s en otro sector más 
pequeño de Tijuana, Baja California (gasto 
medio igual a 1.14 1/s, correspondiente a 123 
casas). La diferencia entre el gasto a cada 60 s 
y el promedio horario es mucho mayor en el 
segundo caso (con el menor número de casas). 
Esta diferencia será aún mayor si el gasto se 
midiera cada segundo, algo que puede ser 
problemático de realizar con los medidores 
de gasto, pero se obtiene directamente con el 
método que se propone en este artículo, permi- 
tiendo introducir en los estudios y proyectos de 
agua potable el concepto de coeficiente de gasto 
máximo instantáneo, algo que por el momento 
sólo se maneja en los estudios y proyectos de 


Variación del gasto 
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Figura 14. Serie medida de gasto cada 60 s (derecha) y su promediado horario (izquierda). Gasto medio = 33.01 1/s; 3 300 casas. 
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Figura 15 Serie medida de gasto cada 60 s (derecha) y su promediado horario (izquierda). Gasto medio = 1.14 1/s; 123 casas. 


ISSN 0187-8336 


Tzatchkov € Alcocer-Yamanaka, Modelación de la variación del consumo de agua potable con métodos estocásticos 


alcantarillado. Se puede esgrimir también el 
concepto de gasto mínimo instantáneo. 

Se llevaron a cabo corridas con el procedi- 
miento propuesto en este artículo, para porcen- 
taje de fugas igual a 18%, tasa de llegada de los 
pulsos de demanda deducida de la variación de 
la demanda de la figura 14, y diferente número 
de casas, desde una sola hasta 100 000 casas, 
y se obtuvieron los coeficientes de demanda 
máxima instantánea y demanda máxima hora- 
ria para cada número de casas. Dado que se 
trata de un proceso estocástico, las corridas 
se realizaron 10 veces (10 realizaciones) para 
cada cantidad de casas, y se obtuvieron los 
promedios de tales corridas. Las figuras 16 y 
17 muestran el coeficiente de demanda máxi- 
ma instantánea y demanda máxima horaria 
obtenidos, comparados con el valor único del 
coeficiente de variación horaria de la deman- 
da que maneja Conagua, que es igual a 1.55 
(MAPAS-Conagua, 2007a). 


Se obtuvieron también los coeficientes de 
demanda mínima instantánea y otros paráme- 
tros. Cabe señalar que los coeficientes obtenidos 
son válidos sólo para los datos (pulsos de de- 
manda y medición del gasto en la fuente) del 
lugar, y pueden ser distintos para otros sitios 
(en particular serán diferentes para la medición 
de la figura 15). Una descripción detallada y 
discusión de los resultados de los coeficientes 
de demanda obtenidos está fuera del alcance 
de este artículo, cuyo objetivo es presentar el 
propio procedimiento de modelación de la de- 
manda, por lo que serán presentados después 
en otra publicación. 


Variación de la demanda en suministro 
intermitente (tandeos) 


El procedimiento arriba explicado es aplicable 
para un servicio continuo de agua potable a 
domicilios sin tinacos o cisternas. Una parte de 
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Figura 16. Coeficiente de demanda máxima instantánea obtenido, comparado con el coeficiente de variación horaria (CVH) 


de la demanda que maneja Conagua. 
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Figura 17. Coeficiente de demanda máxima horaria obtenido, comparado con el coeficiente de variación horaria (CVH) 


que maneja Conagua. 


los domicilios en México cuenta con tinacos o 
cisternas y tinacos, debido al servicio disconti- 
nuo (intermitente) de agua potable. El Instituto 
Mexicano de Tecnología del Agua (IMTA) man- 
tiene un sitio web de Internet de libre acceso 
(http: / /www.pigoo.gob.mx/), creado para 
transparentar la información de una serie de 
indicadores de un gran número de organismos 
operadores de agua potable, alcantarillado y 
saneamiento del país. Uno de los indicadores 
que se maneja en este sitio es el porcentaje de 
tomas domiciliarias con servicio continuo de 
agua potable de 24 horas al día, expresado por 
la relación del número de tomas domiciliarias 
con servicio continuo y el número total de tomas 
en la ciudad. Para el año 2013 participaron con 
información 76 organismos operadores del 
país, que reportaron en promedio un 78% de 
tomas con servicio continuo. Varios organismos 
operadores han realizado esfuerzos importan- 
tes para reducir, o incluso eliminar, el servicio 
de agua por tandeo, mediante programas de 
reducción de pérdidas de agua o aumento de la 
extracción en las fuentes de agua. Se ha notado, 
no obstante, que la mayoría de los usuarios del 
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servicio prefiere mantener las cisternas y tinacos 
en sus casas aun cuando en las tuberías de la red 
se tiene suministro continuo. Algunos aspectos 
hidráulicos del servicio intermitente de agua 
potable, y su modelación, se presentan en Ca- 
brera-Bejar y Tzatchkov (2009), y Cabrera-Bejar 
y Tzatchkov (2012). La variación de los gastos en 
las tuberías con suministro intermitente, o con 
servicio continuo, pero con cisternas y tinacos 
en los domicilios, es muy diferente de aquella en 
tuberías con suministro continuo sin cisternas 
y tinacos. El concepto coeficiente de variación 
horaria de la demanda pierde su significado 
en el caso del suministro intermitente, donde 
la población recibe el agua en horarios preesta- 
blecidos. Dependiendo de la topografía, la red 
de tuberías puede vaciarse cuando concluye 
el horario de servicio y se volverá a llenar al 
reanudar el servicio. En el llenado presentará 
un flujo a tubo parcialmente lleno, cuya mag- 
nitud dependerá de la topografía y desniveles 
existentes, hasta que el agua que entra llene las 
tuberías y expulse el aire. Una vez llenas de 
agua las tuberías, el gasto en las tuberías será 
aproximadamente uniforme. El coeficiente de 
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variación, para fines de revisar si la capacidad 
de la tubería es suficiente, en este caso se define 
por la duración del servicio, sin considerar el 
periodo con flujo a tubo parcialmente lleno. Si el 
gasto medio, en litros por hora, se define como 
el volumen de agua suministrado en un día en 
litros entre 24 horas, y el mismo volumen se su- 
ministra en N horas, el coeficiente de variación 
CV es simplemente: 


CV=24 /N (12) 


Para un servicio de seis horas diarias, por 
ejemplo, se tendría un coeficiente de variación 
de 4, mucho más alto que aquel de un servicio 
continuo, que normalmente se maneja del orden 
de 1.40 a 1.60 (MAPAS-Conagua, 2007a). Es de- 
cir, los gastos y velocidades correspondientes en 
las tuberías que operan con servicio intermitente 
son mucho más altos que aquellos del servicio 
continuo. 

La variación de los gastos en las tuberías 
con suministro continuo, pero con cisternas y 
tinacos en los domicilios, es también diferente 
de aquel en tuberías sin cisternas y tinacos. Es 
de esperar que la variación sea menor en estos 
casos, dado que la demanda de agua en el in- 
terior de los domicilios no llega directamente a 
las tuberías de la red, sino a través de los tinacos 
o cisternas. El problema está bajo investigación 
por los autores. 


Conclusiones 


La variación del caudal en una tubería de agua 
potable tiene carácter estocástico y depende en 
gran medida de la cantidad de casas que abas- 
tece dicha tubería. En este artículo se presenta 
una metodología de simulación numérica, 
basada en el método PRP con inversión de la 
función de la tasa de llegada acumulada, para 
modelar la variación diaria estocástica de la 
demanda instantánea de agua potable, aplica- 
ble para una sola casa o para un grupo de N 
casas, con base en los parámetros estadísticos 
del consumo de agua en las casas individuales, 


nivel de fugas dado, y medición de la variación 
del gasto medida en la fuente de abastecimien- 
to. La metodología propuesta ofrece una base 
racional para la determinación de la variación 
de la demanda, que puede sustituir o combi- 
narse con la manera tradicional de obtener esta 
variación, que se basa solamente en mediciones 
en fuentes de suministro. Los resultados se 
compararon con mediciones de campo en una 
zona aislada de una ciudad mexicana, donde 
el suministro de agua potable es continuo. Para 
esta zona aislada, con la metodología propues- 
ta, se obtuvo que el coeficiente de variación 
horaria de la demanda varía de manera no 
lineal con el número de casas, con valores de 
3.23 para una sola casa hasta 1.64 para 500 o 
más casas, en contraste con el valor constan- 
te de 1.55 que maneja la Comisión Nacional 
del Agua. Se introduce también el concepto 
de coeficiente de variación instantánea de la 
demanda de agua potable, cuyos valores se 
obtuvieron desde 30.58 para una sola casa 
hasta 1.66 para 100 000 o más casas. Por otro 
lado, la variación de los gastos en las tuberías 
con suministro intermitente, o con servicio 
continuo, pero con cisternas y tinacos en los 
domicilios es muy diferente de aquella en 
tuberías con suministro continuo sin cisternas 
y tinacos. El coeficiente de demanda máxima 
horaria es más alto en suministro intermitente 
que en suministro continuo, pero más bajo en 
redes que tienen servicio continuo y cisternas, 
o tinacos en los domicilios. 
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